Anomalous elasticity in nematic and smectic elastomer tubule 
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We study anomalous elasticity in the tubule phases of nematic and smectic elastomer membranes, 
which are flat in one direction and crumpled in another. These phases share the same macroscopic 
symmetry properties including spontaneously-broken in-plane isotropy and hence belong to the same 
universality class. Below an upper critical value D c = 3 of the membranes' intrinsic dimension D, 
thermal fluctuations renormalize the elasticity with respect to elastic displacements along the tubule 
axis so that elastic moduli for compression along the tubule axis and for bending the tubule axis 
become length-scale dependent. This anomalous elasticity belongs to the same universality class as 
that of d-dimensional conventional smectics with D taking on the role of d. For physical tubule, 
D — 2, this anomaly is of power-law type and thus might by easier to detect experimentally than 
the logarithmic anomaly in conventional smectics. 

PACS numbers: 61.30.-v, 61.41.+e, 64.60.F- 



Liquid crystal elastomers [lj are unique materials that 
combine the rubber elasticity of polymer networks with 
the orientational properties and rich phase behavior 
of liquid crystals [2| which includes nematic, smectic- 
A (SmA) and smectic-C (SmC) phases. Elaborate 
crosslinking techniques have been developed to synthe- 
size monodomain samples of nematic Q and smectic [I[ 
elastomers. With these techniques, one can efficiently 
produce from small amounts of material samples in the 
form of thin films or membranes. For example, smectic 
elastomer films have been produced as thin as 75 nm [B| . 
Experiments on such films include measurements of the 
electroclinic effect in flat films [H, 0] and measurements of 
elastic constants of films that have been inflated to spher- 
ical bubbles 0] ■ The potential of liquid crystal elastomer 
membranes for further experimental realizations appears 
promising and calls for a deepening of their theoretical 
understanding. 

Isotropic polymerized membranes have been exten- 
sively studied over the past two decades Q . For example, 
it is well established that a flat phase with long-range 
orientational order in the local membrane normal is fa- 
vored at low temperature over a crumpled phase which 
is entropically preferred at high temperature. More re- 
cently, it has been discovered that permanent in-plane 
anisotropy modifies the phase diagram (jj Hol |: it leads 
to intermediate phases between the usual flat and crum- 
pled phases, so-called tubule phases, which are extended 
in one direction and crumpled in another. Very recently, 
liquid crystal elastomer membranes (see Fi g, p] for car- 




toons) have gained some interest ll|, I12L 1 1 31 ] - in part 
because of their potential to realize anisotropic mem- 
branes experimentally. However, their mesogenic com- 
ponent not only allows for anisotropy, it also allows for 
a spontaneous development thereof. This unique fea- 
ture sets them apart from permanently anisotropic mem- 
branes and provides for a number of interesting phe- 
nomena. Compared to permanently anisotropic mem- 



FIG. 1: Cartoons of elastomer membranes in the (a) SmC- 
flat, (b) nematic-flat, and (c) nematic- or SmC-tubule phase. 
In (a), the thickness of the membrane is vastly exaggerated 
to allow for a depiction of the mesogens and the arrows sym- 
bolize the c-director, i.e., the components of the Frank direc- 
tor perpendicular to the local membrane normal. In (c), the 
arrows symbolize ambiguously the c-director or the director 
depending on whether the membrane is smectic or nematic. 



branes, their phase diagrams are richer, at least if one 
assumes, as we do, idealized crosslinking that avoids 
locking-in permanent anisotropy. The mean-field phase 
diagrams of nematic and smectic (l3| elastomer mem- 
branes each feature five phases, namely isotropic-flat, 
isotropic-crumpled, nematic-flat, nematic-crumpled and 
nematic-tubule for the former, and SmA-flat, SmA- 
crumpled, SmC-flat, SmC-crumpled and SmC-tubule 
for the latter. Because of spontaneous breaking of in- 
plane isotropy, the nematic-flat, nematic-tubule, SmC- 
flat, and SmC-tubule phases exhibit soft elasticity, i.e., 
certain elastic moduli vanish as mandated by the Gold- 
stone theorem, that is qualitatively distinct from the 
elasticity of the flat and tubule phases of permanently 
anisotropic membranes. 

In this note, we study the effects of thermal fluctua- 
tions on the elasticity of the tubule phases of nematic and 
smectic elastomer membranes. These phases share the 
same macroscopic symmetries including spontaneously 
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broken in-plane isotropy and the resulting softness and 
hence belong to the same universality class. Typically, 
fluctuation effects are strong in soft phases because fluc- 
tuations drive elastic nonlinearities, which are often neg- 
ligible in systems without soft elasticity, to qualitatively 
modify the elasticity through a Grinstein-Pelcovits (GP)- 



with the components of the strain tensor u given by 



type renormalization 14] . As a consequence of this 



renormalization, the elasticity becomes anomalous with 
length-scale dependent elastic constants. We explore this 
anomalous elasticity of nematic and smectic elastomer 
tubule by using renormalized field theory . 

Physical membranes are generically two-dimensional 
manifolds embedded in three-dimensional space. In the 
following, to facilitate field theory, we consider gener- 
alizations to D-dimensional manifolds in d-dimensional 
space. For simplicity, we ignore the effects of self- 
avoidence and heterogeneities such as random stresses. 
We employ the framework of Lagrange elasticity the- 
ory [161, [l7j. We label mass points in the unde- 
formed membrane by a reference space vector x = 
(afi, ■ • ■ , %d)- We denote the position in target space of 
the mass point with intrinsic coordinate x by i?(x) = 
(-Ri(x),--- ,i?d(x)), and label the corresponding coor- 
dinates with indices from the middle of the alphabet, 
i,j = l,-",d. To keep our discussion as simple as pos- 
sible, we use orthonormal target space basis vectors ej 
with components e^j = <5y satisfying • ij = Sij and 
choose the reference space basis vectors to form a subset 
of the set {ii} as we can, because the reference space can 
be viewed as a subspace of the target space. 

Physical nematic and smectic tubule phases are charac- 
terized by an equilibrium metric tensor with one positive 
eigenvalue, say £ 2 , and one vanishing eigenvalue [12Ul3l]. 
To facilitate field theory, we consider in the following the 
more general case that there are {D — 1) vanishing eigen- 
values instead [18J. Choosing our basis so that e x = e\ is 
along the eigenvector associated with £ 2 (which is in the 
direction of the nematic or c-director, respectively), we 
can represent the reference conformation of the tubule 
as R° (x) = £ x e x . We label the components of x in the 
reference-space directions perpendicular to e x with in- 
dices from the Greek alphabet, a, (3 = 2, ■ • ■ ,D. Unless 
stated otherwise, the summation convention on repeated 
reference and the target space indices is understood. To 
describe distortions, we use a one-dimensional elastic dis- 
placement field u(x) and a (d— l)-dimcnsional height field 
h(x) which is perpendicular to the tubule's axis e x . With 
this parametrization, the target space coordinate of the 
mass point x after distortion becomes 

i?(x) = [(x + u(x)]e x + h(x). (1) 

The corresponding metric tensor g has the components 



Uxx 
T^xa 



= ±{2£d x u+(d x u) 2 + d x h-d x h}, (3a) 
= \{{C, + d x u)d a u + d x h ■ d a h}, (3b) 
u a /3 = \{d a udpu + d a h ■ dph}. (3c) 

There are various possibilities to set up elastic energy 
densities for liquid crystal elastomer membranes includ- 
ing formulations which account for liquid crystalline de- 
grees of freedom explicitly [l2j, [l3|. For our purposes, the 
most economical one is an effective formulation in terms 
of elastic degrees of freedom only. As a starting point for 
our theory, we use the well established stretching energy 
density for isotropic polymerized membranes [19( which 
we augment with higher-order terms to ensure mechani- 
cal stability under the development of a tubule phase, 



/ = ttrg+ \Btv 2 g + fitig 2 - Ctxgtrg 2 + Etr 2 g 2 



(4) 



where t is a strongly temperature-dependent parame- 
ter, B and /i are Lame coefficients, and C and E are 
higher-order expansion coefficients. In Eq. Q, we omit 
terms beyond fourth order that are irrelevant for our pur- 
poses and we disregard third and fourth order terms such 
as tr 3 g and tr 4 g whose inclusion does not change our 

findings qualitatively. Next, we expand Eq. (j4|) about 
the equilibrium metric tensor describing an undistorted 
tubulus. Dropping inconsequential constant terms we ob- 
tain 



/ = a\u xx + a 2 U aa + biu xx + b 2 u xx u aa + b 3 u 2 aa 
+ h { u lp + 2u lo] - 8Ctr^tr^ 2 + 32( 2 Eu xx trj 
+ 16£tr 2 u 2 . 



(5) 



The elastic constants featured in Eq. ([5]) are conglomer- 
ates of the original clastic constants featured in Eq. (j4|) 
and powers of Q. They are not independent; they obey 
the relations 

a x -oa -C 2 &4 = 0, 6 2 -2fe 3 + 8C 2 C = 0, (6a) 
bi - b 2 + b 3 - 6 4 - 16C 4 £ = . (6b) 

Because u xx describes deviations of g xx from its equi- 
librium value C 2 j the coefficient a± of the term linear in 
u xx must vanish. Exploiting the Ward identities ([6]) and 
setting ai = 0, we can recast the elastic energy density 
as 

/ = 61 wl x + b 2 w xx w aa + b 3 wl a + 64 { - ( 2 u aa + u 2 af3 
+ 2u 2 xa - 2C 2 u xx tr %£ - C 4 u xx tr 2 g 2 } , (7) 

where we have introduced the composite strains 



,9.« 



c 2 



2Un 



9xa — 2u xa , g a p = 2u a/ 3 , (2) 



VJxx — V*xx C tr U , W a a — "ad 



C 2 tru 2 . (8) 
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The &4-term and w aa are of a peculiar structure such 
that multiple cancellations occur once we express / in 
terms of u and h via Eq. ([3]), 

w xx = Cd x u + | (d a u) 2 + irrelevant , (9a) 
w aa = + irrelevant , (9b) 

64-term = — \ 64 ( 2 d a h ■ d a h + irrelevant , (9c) 

where we, anticipating results of power-counting argu- 
ments to be presented below, dropped all terms that turn 
out being irrelevant in the sense of the renormalization 
group (RG). Note that among the various relevant contri- 
butions to the 64-term that cancel, are the contributions 
quadratic in d a u, and thus the tubules are soft with re- 
spect to <9 a it-deformations. In order to proceed towards 
the desired field theoretic Hamiltonian, we insert Eq. (|9|) 
into the elastic energy density (J7J) and we rescale x — > Qx. 
Then we add bending terms as mandated by mechanical 
stability and finally integrate over the reference space 
coordinate x to switch from the elastic energy density. 
These steps result in 



n 

T 



1 

2T 



dxl B 



1 



d x u + -(d a u) 



+ B h d a h ■ d a h + K u (d 2 a u) 2 + K h d 2 x h -dlh\, (10) 



where B u = 2bi, Bh = — (£ 2 /2)&4 and where K u and 
Kh are bending moduli. When expressed in terms of 
the original elastic constants and (, Bh reads Bh = 
4£ 4 [C— ( 2 E— fi/C 2 }- For mechanical stability, Bh must be 
positive implying that C has to satisfy C > Q 2 E + /i/C 2 - 

At this point it is worthwhile to discuss the rotational 
symmetries of the Hamiltonian (|10[) briefly. As is the case 
for any elastic medium, the elastic energy describing our 
tubule must be invariant under global rotations in target 
space. Moreover, because we are interested in tubules 
that emerge from a phase with no in-plane anisotropy 
via spontaneous symmetry breaking, our elastic energy 
must also be invariant under global rotations in reference 
space. However, because we neglect irrelevant terms, Ji 
can satisfy these symmetries at best for small rotation 
angles, and indeed it does, as can be checked straightfor- 
wardly, e.g., by rotating the reference and target space 
bases using appropriate rotation matrixes. 

Next, we turn to the aforementioned power-counting 
analysis to assess the relevance of the various terms in 
the sense of the RG. To have some guidance in the anal- 
ysis, we first consider the mean square fluctuations of 
the displacement and height fields in real space in the 
harmonic approximation. The Hamiltonian (|10|) implies 
that the Gaussian propagators of these fields are given in 
momentum or wave-vector space by 

^fa) = 5 I2 T is „4 > G h,ij(q) = 



B u q 2 x + K u q\ 



B h q\ + K h q x ' 
(11) 



with q\ = q 2 a and where it is understood that % and j run 
over the subspace of the target space that is perpendicu- 
lar to the tubule's axis e x . Calculating the Gaussian fluc- 
tuations of the displacement field in real space by Fourier- 
transforming G u (q), we find (u(x)w(x)) ~ L 3 ^ D , where 
L±_ is the length of the tubule in any of its directions 
perpendicular e x were it to be flattened out. Thus, in 
the harmonic approximation, the u-fluctuations diverge 
in the infra-red for D < 3. For the height field, on the 
other hand, we have (/i(x) • h(x)) ~ L 5 / 2 D , which di- 
verges in the infra-red for D < 5/2. Thus, if we decrease 
D from a high value where the mean field approximation 
is exact to lower dimensions, infra-red divergences start 
to occur at D = 3, which signals that D c = 3 is the upper 
critical dimension. This observation in conjunction with 
a inspection of G u also signals that in our power-counting 
analysis we should count each power of q x as two powers 
of q a . Note in comparison, that tubules in permanently 
anisotropic membranes do not feature the soft elastic- 
ity with respect to 9 Q u-deformations, and that their it- 
propagator is thus qualitatively different from that given 
in Eq. (fTT]l in that it carries a q\ instead of the q\_ 
As a consequence, the real-space Gaussian fluctuations 
in the latter become divergent in the infra-red in D < 2 
rather than in D < 3 leading to an upper critical dimen- 
sion D c = 5/2 rather than D c = 3. As far as the lower 
critical dimension Di c is concerned, i.e. the dimension be- 
low which the membrane is inevitably crumpled, inspec- 



tion of harmonic fluctuations (d x h(x) ■d x h(x)) ~ L" ± 



-3/2-D 

1 

of the tubule normal along e x indicates that Di c = 3/2. 

Now, we determine which terms are relevant in the 
sense of the RG by rescaling our coordinates such that 



the resulting coordinates are dimensionless: x 



H x 



and x a — > [i~ x a , where \i is an inverse length scale that 
must not be confused with the Lame coefficient encoun- 
tered earlier. Under this rescaling, the Hamiltonian (fT0|) 
remains invariant in form provided that h — * /ih, T 
fi 3 ^ D T and K h -> ^ 6 K h . This means that h, T and K h 
have the naive or engineering dimensions 1, e = 3-D 
and —6, respectively. The field u and the remaining pa- 
rameters in Eq. (|10|) have a vanishing naive dimension. 
Above D — 3, the naive dimension of the temperature 
T is negative implying that T is irrelevant above D = 3. 
This signals that D c = 3 is the upper critical dimension 
in accord with what we have seen above. The Hamil- 
tonian (|10[) contains all relevant terms. In addition, it 
contains the if/j-term although Kh has a negative naive 
dimension. Dropping the A'/i-term, would make the h- 
propagator independent of q x which is un-physical, i.e., 
Kh is a dangerous irrelevant coupling constant that must 
be kept. 

The Hamiltonian (fT0|) has the remarkable feature that 
there are no relevant contributions that couple the dis- 
placement and the height field. Thus, it decomposes into 
a part that depends only on u and a part that depends 
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only h. The ft,-part is purely harmonic. Therefore, there 
is no anomalous elasticity with respect to h in the vicin- 
ity of D c . The u-part is equivalent to the Landau-Peierls 
Hamiltonian as studied by GP with d replaced by D. 
Consequently, there is anomalous elasticity with respect 
to u and this anomalous elasticity belongs to the GP uni- 
versality class with D taking on the role of d. For bulk 
smectics as studied by GP, the physical and the upper 
critical dimensions coincide, and hence the compression 
and bending moduli B u and K u depend logarithmically 
on wave-vectors. For our tubule, however, the physical 
case is D = 2 < D c , and hence B u and K u have power- 
law behavior, 

B u (ci) = q2 B S B (q x /\q ± n, (12a) 
K u (q) = qrS K (q x /\q ± n, (12b) 

with critical exponents given to one- loop order by i]b = 
|e, tik = — \s and z — 2 — |s (e = 1 corresponds to 
the physical case). Sb and Sk are scaling functions with 
the asymptotic properties Sb(u) ~ Sb(v) ~ const for 
y — * oo and Sb(v) ~ and Sk(v) ~ V ~ VK for y — > 0. 

The signs of t\b and t\k imply that B u and K u , respec- 
tively, vanish and diverge at long length-scales. Thus, 
just as the interplay of thermal fluctuations and elastic 
nonlinearities stabilizes the flat phase in isotropic poly- 
merized membranes [20| . it here stabilizes the tubule 
phases through infinitely enhancing the bending modulus 
K u at long length-scales. 

An important question is whether our results remain 
valid for D < 5/2. On one hand, Gaussian theory in- 
dicates that fluctuations of h become divergent in these 
dimensions. On the other hand, it is entirely possible 
that the system is deep enough in the interacting regime 
for D w 5/2 that the Gaussian fixed point of the RG has 
lost its significance and that, hence, the true lower critical 
dimension for h may be considerably smaller than indi- 
cated by Gaussian theory. A similar problem arises, e.g., 
in the dynamics of ferromagnets 21|, where the upper 
critical dimension is 6, and, though the upper critical 
dimension for corresponding statics is 4, it is generally 
believed that the dynamical field theory remains useful 
down to the physical dimension 3. To settle this ques- 
tion definitely for the tubule, one has to devise a RG 
study that reliably treats both u and h below D = 5/2, 
which we leave for future work. We hope that our work 
stimulates further experimental and theoretical interest 
in liquid crystal elastomer membranes and their tubule 
phases. It would be interesting to see measurements of 
their elastic moduli, e.g., by light scattering, sound or 
stress-strain experiments. Because the dependence of B u 
and K u on wave- vectors is of power-law type rather than 
logarithmic, their anomalous elasticity might be easier to 
detect than that of conventional smectics. 
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